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ABSTRACT 

It has been known for more than 30 years that star formation in giant molecular clouds (GMCs) 
is slow, in the sense that only ~ 1% of the gas forms stars every free-fall time. This result is 
entirely independent of any particular model of molecular cloud lifetime or evolution. Here we survey 
observational data on higher density objects in the interstellar medium, including infrared dark clouds 
and dense molecular clumps, to determine if these objects form stars slowly like GMCs, or rapidly, 
converting a significant fraction of their mass into stars in one free-fall time. We find no evidence for a 
transition from slow to rapid star formation in structures covering three orders of magnitude in density. 
This has important implications for models of star formation, since competing models make differing 
predictions for the characteristic density at which star formation should transition from slow to rapid. 
The data are inconsistent with models that predict that star clusters form rapidly and in free-fall 
collapse. Magnetic- and turbulence-regulated star formation models can reproduce the observations 
qualitatively, and the turbulence-regulated star formation model of Krumholz & McKce quantitatively 
reproduces the infrared-HCN luminosity correlation recently reported by Gao & Solomon. Slow star 
formation also implies that the process of star cluster formation cannot be one of global collapse, 
but must instead proceed over many free-fall times. This suggests that turbulence in star- forming 
clumps must be driven, and that the competitive accretion mechanism does not operate in typical 
cluster-forming molecular clumps. 
Subject headings: ISM: clouds — stars: formation 



1. INTRODUCTION 

iZuckerman fc Evand l)1974|) first pointed out that star 
formation in giant molecular clouds (GMCs) happens 
surprisingly slowly. Comparing the mass of GMCs in the 
Galaxy with the total Galactic star formation rate im- 
plies that no more than ~ 1% of the gas can form stars 
for each cloud free-fall time. This result is sufficiently 
surprising that numerous theories have been proposed 
to explain it, ra nging from the ide a that strong mag- 
netic fields (e.g. Allen & Shu '200Cf) or turbulence (e.g. 
^rumholz & McKce 2005) within clouds inhibit star for- 
mation to the idea that galactic-sca le gravitationa l insta- 
bility regulates star formation (e.g. lLi et al.l2005aj) to the 
idea that GM Cs are, contrary to most observational esti- 
mates to date fBlitz et al."2005") , actually gravitationally 
unbound (e.g. Clark & Bonnell 20042. 

FoUowing lKrumholz fc McKed l)2005j) , we define the di- 
mensionless star formation rate per free-fall time SFRff 
as the fraction of an object's mass that it converts into 
stars per free-fall time at the mean density of that ob- 
ject. The Zuckcrman & Evans (1974) argument shows 
that SFRff « 0.01 for GMCs. This provides a powerful 
constraint on models of GM Cs. For e xample, i t rule s 
out the early GMC model ofl^dreich fc KwanI ljl974j) . 
a simplified version of which is that GMCs are spheres of 
gas of density p that are in free-fall collapse. The clouds 
reach a singularity in a time iff(/5) = [37r/(32G'/9)]^/^, at 
which point their mass is converted into stars. For a pop- 
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ulation of such clouds SFRff = 1, w hich is inconsistent 
with the ^Zuck erman fc EvansI l)1974l) result. 

An important observational question, and a crucial 
test for theories of how star formation is regulated, is 
to what densities and length scales SFRg remains much 
smaller than unity. In other words, is there a d ensity 
at which something like the iGoldreich fc Kwai] l)1974f) 
free-fall collapse model becomes reasonable? As an ex- 
ample, consider observing a star-forming region using a 
molecular tracer sensitive to gas at densities of nn ^ 10^^ 
cm~^, where ur is the number density of hydrogen nu- 
clei, to estimate the total mass of gas at such high den- 
sities. This is larger than the m ean density of "cores" 
seen both observationally (e.g. iBarranco fc GoodmanI 
Il998^ and in simulations of star form ation regulated 
by turbulence (e.g. Jappsen 'et alJl2005|l . and is roughly 
the density at which models of magnetically-regulated 
star formation predict that gas will completely decou- 
ple from the magnetic field and enter free-fall collapse 
fe.g. lDesch fc Mouschovi a^2001'). At such high densities 
protostellar outflows can proba bly stop at most ha lf the 
gas from reaching a star ((Matzner fc McKe3 l2000'^ . and 
the thermal pressure in gas at 10"'^^ cm~'^ is considerably 
higher than the typical ram pressure of the turbulence in 
GMCs, so gas at such high densities is largely impervious 
to external perturbations. 

Thus, regardless of the model of star formation one 
adopts, one would expect that almost all of the gas at 
densities ^ 10^^ cm~'^ is part of gravitationally-bound, 
collapsing objects that have largely decoupled from the 
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background turbulent flow. In the absence of effective 
internal support or external disturbance, order unity of 
the gas at such high densities is likely to be incorporated 
into a star within one free-fall time. For this reason, 
essentially all models of star formation predict that the 
total mass of gas at densities ^ 10^^ cm~^, divided by 
the free-fall time of this gas, should yield a value compa- 
rable to the total star formation rate in the region over 
which the mass is measured. Instead of SFRff 0.01 
as for GMCs, i.e. slow star formation, one would obtain 
SFRff ~ 1, i.e. rapid star formation. 

However, different models make different predictions 
about the shape of the curve of SFRff versus density in 
between ~ 1% at the characteristic density of GMCs, 
~ 100 cm^^, and ^ 1 at a density ^ 10^^ cm"'^. These 
different predictions correspond to different models of the 
physical scale at which gas both decouples from the back- 
ground flow and ceases to be supported by internal feed- 
back mechanisms, and thus transitions from slow to rapid 
star formation. At one extreme, magnetic regulatio n 
models such as those of iDesch fc MouschoviasI ()2001() . 
neglecting for the moment turbulent enhance ment of the 
ambipolar diffusion rate (e.g. "lleits ch et alil2004D . pre- 
dict that star formation only becomes rapid once gas 
decouples from the magnetic field, a process that does 
not even begin un til densities of ^ 3 x 10^° cm~^. At 
the other extreme, iBonnell et alJ l|2fl03D argue that star 
clusters form from gas clumps at densities of ~ 5 x 10'* 
cm~^ that undergo a free-fall collapse in which at least 
~30% of their mass is converted into stars l|KrouDa et al.l 
120011) . so star formation should be rapid at this density 
or higher. In this model, the decoupling scale corre- 
sponds to the transition from globally unbound struc- 
tures (GMCs) to globally bound structures (pro t oclus- 
ters). Thus, extending the iZuckerman fc Evanij l|1974f) 
calculation of star formation rate divided by free-fall 
time to higher densities, in hopes of identifying a scale at 
which there is a transition from slow to rapid star forma- 
tion, provides a means of distinguishing between models 
of how star formation is regulated. 

It is important at this point to differentiate the con- 
cepts of the rate and efficiency of star formation, and 
the lifetime of star-forming clouds. Unfortunately these 
terms are often confused in the literature, and there 
are no standard definitions, so we describe here the 
definitions we use in this paper. The star formation 
rate M» is the easiest to define, since it is simply the 
instantaneous conversion rate of gas into stars within 
some volume. If we pick a density threshold p, we can 
then define the dimensionless star formation rate per 
free-fall time for the gas above that density threshold, 
SFRff = M,/[M(> p)/ts{p)], where M(> p) is the mass 
of material inside the volume of density p or greater. 
In contrast, the lifetime ici(p) of star- forming cloud is 
somewhat more ambiguous. We take it to mean the to- 
tal duration during which a cloud is visible in a tracer 
sensitive to densities of p or more. Note that our defini- 
tion neglects the complication that something that starts 
as a single cloud may in the course of its evolution break 
up into multiple pieces, so the visible lifetime and the 
dynamical lifetime may be different. Finally, by the effi- 
ciency e{p) we mean the fraction of gas mass M(> p) that 
is converted into stars by a cloud over its lifetime tci(p)- 



Again, we neglect the complication that clouds are not 
closed boxes, so M(> p) is likely to be time-dependent 
due to continuing accretion or mass loss. Ideally e{p) 
should be computed using a Lagrangian definition of the 
cloud mass, i.e. all the mass that reached density p or 
more at some point. However, this is generally not a 
direct observable. 

Roughly speaking, the rate, efficiency, and lifetime are 
related by M, ~ e{p)M{> p) / tdoudip) ■ Consequently, 
one can define a rough time-averaged value for SFRff 
in a cloud (SFRff(/r?)) « €{p)/[tci{p)/tff{p)], i.e. the mass 
fraction converted into stars divided by the cloud lifetime 
in free-fall times. One could also describe (SFRff) as 
being the efficiency per free-fall time, since it measures 
a fraction of mass converted into stars. However, we 
refer to it as a rate because it is measured in amount 
per unit time, whereas efficiency in the literature most 
commonly refers to the total fraction of mass converted 
into stars, not the amount per unit time. Moreover, since 
the time-averaged definition is ambiguous anyway due to 
uncertainties in exactly what is meant by the efficiency 
and the lifetime, we will generally use the instantaneous 
value of SFRff, which is well defined and, as we show, 
directly observable. 

The reason for making all these definitions explicit 
it to point out that observational constraints on SFRff 
by themselves do not directly constraint ^ci or e, and 
vice versa. Numerous authors have used various ob- 
servational t echniques to estimate tr] in clusters and 
GMCs (e.g . Elmegreen 2000; Hartmann et al. "200^ 
Tassis fc Mouschov ias 2004; Mouschovia s et al. 2005t 
Tan et alJ l2006HBa llesteros-Pa.redes fc H artmannI feOOa 
Blitz et al. l2005j) . However, with the exception of S 13.51 
we will not discuss the question of cloud lifetime at all 
in this paper. When we refer to rapid versus slow star 
formation, what we mean is SFRff ~ 1 or ^ 1, not is 
tci ^ tff or :s> tff. These questions a re conceptually dis- 
tinct. As an example, note that in the lGoldreich fc KwanI 
ljl974[) picture of GMCs as collapsing spheres, the evo- 
lution time scale is always the cloud free-fall time, td 
tg. However, the model still gives SFRff — 1, and we 
would ther efore describe it as a rapid star f orma tion. 
In contrast. IBallesteros-Paredes fc HartmannI l|200(a) ar- 
gue that star formation lasts only one crossing time of a 
molecular cloud, so again td ^ iff, but that during this 
time only 1% of the mass turns into stars. We would 
describe this as slow star formation, since SFRff ^ 0.01, 
eve n though the cloud e volution time is simila r to that in 
the IGoldreich fc KwanI model. Th e Zuckerm an fc EvansI 
argument rules out the IGoldreich fc KwanI free-fall col- 
lapse model for GMCs, but i s consistent with the 
IBallesteros-Paredes fc HartmannI model. Thus, we em- 
phasize that in this paper we remain completely agnostic 
on the question of molecular cloud lifetime. 

In this paper we consider star formation in several 
classes of object. Infrared dark clouds (IRDCs) are 
regions of high extinction seen in a.bsorption against 
the Galactic infrared background l)Egan et al.l 119981: 
iCarev et al.ll2000t iSimon et al.ll200^ . IRDCs are clearly 
associated with star formation, and in at least some cases 
IRDCs have ma ssive protostars embedded within them 
iRathborne et al. 2005). S evera l authors have suggested 
(e.g. IMente~t al.. .20051 ITmI l2005l IRathborne eTal] 
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1200^ that IRDCs are the progenitors of star clusters. 
Within IRDCs, at still higher densities, are dense molec- 
ular clumps. These objects may be observed in a vari- 
ety of molecular transitions with high critical d ensities, 
and we co nsider tw o here: HCN(l-O) (Gao & Solomon! 
|2hn4a.b„: Wu et al 120051 and CS(5-4) (Plume et a l. 1997i: 
iShirlev et al.ll2003j) . Molecular clumps seen in these two 
transitions are often associated with water masers and 
other signs of massive, clustered star formation. Our 
goal is to determine SFRff for each of these increasingly 
dense gas tracers. We also determine this quantity for 
the Orion Nebula Cluster using a completely different 
method, which provides an independent check on our es- 
timates. 

The remainder of this paper proceeds as follows: in §|5] 
we use a variety of observations to derive SFRff for our 
objects and construct a plot of SFRff versus character- 
istic density to search for signs of a transition from slow 
to fast star formation. In § Ol wc compare our results to 
theoretical models for the regulation of star formation, 
and also point out some implications for the process of 
star cluster formation. Finally, in §0]we summarize our 
conclusions, and suggest directions for future work. 

2. ESTIMATES OF SFRff 

Let Af, be the total star formation rate in a galaxy, 
and consider star formation occuring in objects of class 
X. Let fx be the fraction of galactic star formation that 
occurs in these objects, Mx be their total mass in the 
galaxy, and tg-x be their typical free-fall time. Then 



In this section we use observations to estimate all the 
factors on the right hand side, and thus determine the 
dimensionless star formation rate SFRff_x in several 
classes of object. We emphasize that, as we vary the 
typical density of the objects we consider, both tg-x and 
Mx will vary - as we go to denser and denser clouds, the 
free-fall time will decline, pushing SFRff_x downward, 
but the amount of mass in objects of that density will 
also decrease, pushing SFRff back up. Ultimately, we 
wish to determine which of these two effects dominate, 
or whether they roughly balance. 

In addition to determining a best value for SFRff _x, 
we also attempt to make a rough estimate our uncertain- 
ties by estimating the uncertainty on each parameter in 
equation We do this from data where possible and 
from more general considerations where not, and we add 
the uncertainty factors in quadrature. This method of 
combining errors is formally appropriate only if the er- 
rors are independent and Gaussian-distributed, neither 
of which are likely to be completely true, but the proce- 
dure does give us a sense of by how much our estimates 
could vary. 

2.1. Star Formation Rates 

The unknown on the right hand side of that has 
been studied most heavily is the star formation rate. 
Numerous authors have discussed methods of inferring 
the star formation rate from various observ ables (e.g. 
IKenni cutt 1998a; l erlesias-Paramo et al. 2006, and refer- 
ences therein), so here we only summarize aspects of this 
discussion that we will apply directly in what follows. 



2.1.1. The Milky Way 
In the Milky Way, iMcKee fc Williamsl ifTool find 
Af, w 3 Mq yr~^ based on catalogs of Galactic HII re- 
gions. Star formation is distributed in an exponential 
disk with scale radius — 3.3 kpc and sharp cutoffs at 
3 kpc and 11 kpc in Galactocentric radius. For this distri- 
bution, approximately 2/3 of Galactic star formation oc- 
curs within the solar circle. Th e dominant uncert ainty in 
this estimate, roughly 0.3 dex l|Kennicut ti 1 1 9 9 8 aj) . is the 
shape of the stellar initial mass function (IMF), since HII 
regions only trace the massive stellar population and one 
must extrapolate to estimate the total mass. Estimates 
of the star formation rate based on chemical modelling 
(jPrantzos fc Auber3ll995|) also give star formation rates 
consistent with 3 Mq yr~^ to within a factor of 2. 

2.1.2. Extragalactic Far-Infrared Observations 

For external galaxies, a commonly-used tracer of star 
formation is far- infrared (FIR) light. While FIR is 
usually not the preferred tracer of star formation in 
normal spiral galaxies like the Milky Way, it is de- 
tectable over a very wide range of sources, from ultralu- 
minous infrared galaxies (ULIRGs) to individual cluster- 
forming clumps in the Milky Way, making it possible to 
look for correlations over a wide range of l uminosities 
(e.g. 'Gao fc Solomon"2004b': lWu et al.feoO.'iD . iKennicutii 
|l998a) finds that in optically-thick starbursts the IR lu- 
minosity and star-formation rate are related by 

AU/ {Mq yr-1) = 2 x IQ-^^Lm/LQ). (2) 

The dominant uncertainty is the age of the stellar pop- 
ulation, but as we discuss in more detail below, for ex- 
tragalactic sources this is only ^ 30%. The uncertainty 
is considerably larger in normal spiral galaxies, where 
old stellar populations contribute a significant luminos- 
ity, and where only a fraction of the light is reprocessed 
into infrared. Based on comparisons of multiple tracers 
of the star formation rate in a range of galaxy types, 
llglesias-Paramo et alJ ()2006j) find that ^ probably un- 
derestimates the star formation rate in normal spirals. 
We instead adopt 

M,/ {Mq yr-1) = 4 x 1O-i°(Lir/L0) (3) 

in normal spirals, which is approximately consistent with 
the mean in the sample of llglesias-Paramo et al.l l|2006fl 

for low star-formation rate (A^, ^ 10 Mq yr~^) systems. 
There is, however, a much larger scatter in this relation 
than in the corresponding relation for starburst systems. 

2.1.3. Galactic Far-Infrared Observations 

Several authors (e.g. "Plume et al. 1997"; IMueller et alJ 
1^02; Shirlcv ct al. 2003; Wu ct al. 2005) use infrared lu- 
minosities to estimate the star formation rates within 
individual cluster-forming gas clumps in the Milky Way. 
These objects have very high surface densities, so essen- 
tially all the light is reprocessed into IR and the uncer- 
tainty for spiral galaxies does not apply. However, we 
will not directly compare to these data because there 
is a much larger uncertainty arising from the age of the 
stellar population. In systems 3 Myr old no stars have 
disappeared through supernovae, so the massive stellar 
population cannot yet have reached equilibrium between 
formation and destruction. 
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Fig. 1. — Luminosity to star formation rate ratio (upper 
panel) and luminosity to stellar mass ratio (lower panel) versus 
time for a stellar population forming at a constant rate. We 
show the f ollowing computations: using starburst99 version 5.0 
|Lcithcrcr ct al. 1999; Vazaucz & Lcithcrcr 2005) with the default 
parameters (dashed lines); using starburst99 with an IMF that 
has a slope of —2.35 from 0.1 Mq to 120 A/q {dotted lines); and 
the median (thick solid lines) and Icr upper and lower limits (thin 
solid lines) for 400 clusters simulations computed wit h the m odel 
of .Tan fc McKcc I.2002.V Further details are given in S 12.1.31 



To study the magnitude of the uncertainty that 
this effect induces, w e use starburst99 version 5.0 
l|Leitherer et al.l Il999t iVazauez fc Leithered l200fJ) to 
compute the evolution of Hght to star formation rate and 
hght to mass ratios for systems with constant star for- 
mation rates. Figure ^ sho ws the results. The dashed 
lines show models using the iKrouDal ll2 nnl IMF while 
the dotted lines use an IMF with the ISalneteri ^95^) 
slope of a = -2.35 from 0.1 Mq to 120 Mq] for all other 
starburst99 parameters, we use the defaults. As the plot 
shows, at ages of a few Myr or more, the luminosity per 
unit star formation rate is quite insensitive to both age 
and IMF, so luminosity is a good indicator of star for- 
mation rate. However, at younger ages the luminosity 
in the starburst99 calculation traces total stellar mass 
more closely than star formation rate, so the light is not 
a good indicator of star formation rate. 

One might think it possible to break this degeneracy 
by independently determining the age of the stellar pop- 
ulation, and then using the luminosity to infer the stel- 
lar mass and thus the star formation rate. However, a 
more detailed treatment of very young systems than star- 
burst99 provides shows that infrared observations alone 
are not sufficient to constrain the stellar mass in systems 
^ 1 Myr old. In such young systems, a significant frac- 
tion of the stars may not yet have contracted to the main 
sequence, in which case they will be more luminous than 
equal mass stars on the main sequence would be. Accre- 
tion luminosity may further enhance the radiative output 
above what would be found for a non-accreting popula- 
tion of the same mass. There are also factors that reduce 
the luminous output compared to an older population of 



the sam e mass. Massiv e stars require '--^ 0.1 Myr to as- 
semble ()McKee fc Tani r2002. 2003), so systems younger 
than this will be missing their contribution to the light. 
In systems of 1000 Mq or less, poor sampling of the 
massive tail of the stellar IMF produces a large scatter, 
and causes the median system to be less luminous per 
unit mass or per unit star formation rate than a larger 
population would be. 

We explore how these effects ch ange the light output o f 
a young cluster using the model of lTan fc McKel l|200^ . 
In Figure ^ the solid lines show the hght output from 
a simulated cluster with a final stellar mass of 1000 Mq 
forming at a constant star formation rate of 1.08 x 10"'^ 
Mq yr~^. This model uses the same IMF we use for the 
Salpeter-slope starburst99 calculation (the dotted line in 
the figure) , but it includes accretion luminosity, pre-main 
sequence evolution, finite star formation times, and dis- 
crete sampling of the IMF. It does not include any post- 
main sequence evolution. The thick central line shows 
the median of 400 runs with different samplings of the 
IMF, and the thin lines above and below it show the 
tracks that bound 68% of the runs. As the plots show, 
the median cluster ^ 1 Myr in age will be much less lumi- 
nous per unit star formation rate than a galactic stellar 
population that is ^ 10 Myr old and contains enough 
stars to fully sample the IMF, but the spread can be an 
order of magnitude for a 1000 Mq cluster. As the cluster 
mass increases the effects of discrete sampling decrease, 
causing the luminosity spread to decrease and the median 
to rise. Even for clusters massive enough to sample the 
full IMF, though, neither the light to mass ratio nor the 
light to star formation rate ratio stay constant at ages 
^ 1 Myr, so light is a poor indicator of either stellar mass 
or star formation rate. 

While we cannot use infrared luminosity to study 
the star formation rate in Galactic star-forming gas 
clumps directly, there is an observational correlation be- 
tween infrared and molecular luminosity for such objects 
ifWuet al. 2005 ) fr om which we can learn a great deal. 
We discuss how to interpret this correlation in light of 
our results in § 13.51 

Finally, note that there are few star clusters where one 
can estimate the star formation rate directly by plac- 
ing a large number of stars on the HR diagram and 
using pre-main sequence tr acks to estimate the c luster 
mass and age spread (e.g. iPalla fc Stahled 119991 120001: 
Fluff & Stahlcr 200^. Since this requires luminosity and 
temperature determinations for many stars, it is posssi- 
ble only in systems without too much extinction, which 
limits this technique to low density regions (e.g. Taurus) 
or somewhat older regions where most of the initial gas is 
gone (e.g. the Orion Nebul a Cl uster, ONC). We discuss 
the ONC i n more detail in 12.51 and also refer readers to 
iTan et all |2006) for a detailed discussion of other tech- 
niques by which one may estimate the formation times 
of star clusters. 

2.2. Star Formation in Infrared Dark Clouds 

As discussed above, IRDCs are likely to be the pro- 
genitors of star clusters. Simon ct al. ( 2006) catalog the 
infrared dark clouds in the inner Galaxy. This enables 
us to compute the star formation rate per free fall time 
in IRDCs by comparing to the total star formation rate 
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in the inner Galaxy, roughly 2 Mq yr^ with a factor 
of 2 uncertainty. We should note that our calculation 
of the star formation rate in IRDCs is entirely different 
from the one given Rathborne et al. (2006). Their esti- 
mate assumes an unknown time scale and efficiency of 
star formation. The calculation we give here does not 
rely on any such assumptions. 

First we must estimat e the fraction of star formation 
that occurs in IRDCs. iLada fc Lada (2003) find that 
70 — 90% of Galactic star formation occurs in clusters. 
If a cloud were visible during the entire star formation 
process as an IRDC, this would imply /irdc ~ 0.8. 
However, clouds that have too many embedded proto- 
stars will not be infrared dark, and will therefore cease 
to be visible as IRDCs. This probably does not occur 
until most of the stars have formed, though, so we adopt 
/irdc — 0.8 as a reasonable guess, while acknowledging 
that it could be a bit smaller. Being slighly more c on- 
servative about our uncertainty than ILada &: LadaL we 
adopt /irdc = 0.6 — 1 as the plausible range of possibil- 
ities, i.e. anywhere from 60% to 100% of star formation 
inside the solar circle takes place in IRDCs, correspond- 
ing to a factor of 1.25 uncertainty on our estimate of 
/irdc = 0.8. Note that /irdc — 1 gives the largest pos- 
sible value of SFRff for IRDCs, so any uncertainties that 
move /irdc out of our plausible range can only make 
S FRff smaller, not l arger 

iRathborne et all lj2006|) find that the total mass of 
IRDCs in the inner Galaxy is Mjrdc ~ 10* Mq based 
on the measured properties of a sub-sample of IRDCs 
observed in 1.2 mm continuum emission, and an esti- 
mate o^thB_detecti£ii efficiency for the MSX IRDC sur- 
vey l)Simon et alJl200'6|l . This mass estimate is probably 
uncertain by a factor of several, becaus e it is unclear how 
representative the clouds surveyed by ' Rathborne et al.l 
((2006j) are of the entire IRDC population. The sam- 
ple consists of the darkest clouds (darkness measures a 
combination of column density and degree of background 
illumination) from a sample with known kinematic dis- 
tances. This selection introduces an unknown bias in the 
mass estimate, so we consider the IRDC mass estimate 
to be uncertain by factors of a few. 

The free-faU time in IRDCs is tg = [37r/(32Gp)]i/2, 
where p is the mean densit y . For the 38 IRDCs in the 
sample of IRathborne et alJ lj2006|) , we define an effec- 
tive radius r = {A/ttY^^D, where A is the angular area 
within the 2a detection thr e shold of the cloud taken from 
the catalog of lSimon et al.l lj2006(l . and D is the distance 
estimate taken from Rathborne et al. We take the mean 
density to be p = 3M/{ATTr^). This procedure is fairly 
uncertain, since the location of the 2a contour depends 
on the background emission, and the morphology is fila- 
mentary rather than round for a significant minority of 
clouds. Nonetheless, we can make a rough estimate for 
p and tff, and check for any systematic variations with 
IRDC size. We plot the derived free-fall times and den- 
sities in Figure[21 As the plot shows, there is a spread of 
a factor of 3, but no clear systematic trend. We adopt 
as our characteristic free-fall time the mass- weighted har- 
monic mean 

te-iR,DC = ^ = 0.9 Myr, (4) 
corresponding to a number density of hydrogen nuclei 
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Fig. 2. — Free-fall time and density versus mass for the 38 
IRDCs in the sample of lRathborne et al.l 12009) . 



uh = 2 X 10^ cm The range of free-fall times covering 
the central 67% of the IRDCs is 0.45 - 3.0 Myr, so the 
uncertainty is roughly a factor of 2.6. Similarly, the cen- 
tral 67% of IRDC densities covers a range of 210 — 9200 
cm~^, an uncertainty of a factor of 6.6. 

We now have all the necessary figures to plug into equa- 
tion JQ: the fraction of star formation inside the solar 
circle taking place in IRDCs is /irdc = 0.8, with a 25% 
uncertainty; the total star formation rate inside the solar 
circle is 2 Mq yr~^, with a factor of 2 uncertainty; the 
free-fall time in IRDCs is 0.9 Myr, with a factor of 4 un- 
certainty; the total mass of IRDCs in the inner Galaxy 
is 10* Mq, with a hard-to-quantify uncertainty, which 
we take to be a factor of 3. Combining these estimates 
gives a fiducial value SFRff_iRDc ~ 0.014. Combining 
the errors as described in § |21 gives an uncertainty of a 
factor of 4.6. 

2.3. Star Formation in Gas Traced by HON 

The molecular transition HCN fl-O) has a critic al den- 
sity of = 6 X lO^/r cm^^ fc.f. lGao fc Solomonll2004aL 
who give the critical density in terms of number den- 
sity of hydrogen molecules rather than hydrogen nuclei), 
where r is the line-center optical depth in the escape 
probability approximation, and is therefore an excellent 
tracer of the dense molecular r egions as sociated with star 
formation. iGao fc SolomonI l)20n4alH l observe a large 
sample of normal spiral galaxies and luminous and ul- 
traluminous infrared galaxies in HCN(l-O) and show that 
the HCN luminosity correlates well with the IR luminos- 
ity, following 

911 ± 227 ^ , (5) 

Lq K km s pc 2 

where the figure quoted is the mean plus or minus the 
standard deviation. Galactic cores with luminosities 
above Lq observed in HCN by IWu et abl (|200l 

show a similar correlation, with almost exactly the same 
mean, although the median is lower than the median of 
the extragalactic sample by a factor of 2.5. Since in- 
frared luminosity is a tracer of star formation rate, and 
HCN(l-O) is a tracer of molecular gas, this relation is a 
direct measure of the star formation rate per unit mass 
in gas of densities traced by HCN, and thus of SFRff in 
that gas. We use the difference in medians between the 
Galactic and extragalactic samples as an estimate of the 
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uncertainty factor in the IR-HCN correlation, although 
we note that this is almost certainly an overestimate for 
the extragalactic sample because that data set is much 
more tightly correlated than the Galactic sample. 

To derive the total gas mass that corre sponds to a given 
luminosity Lrcn in the HCN(l-O) line. lWu et all (|2QQ^ 
observe in HCN(l-O) a sample of star-forming clumps 
with known virial masses determined from optically-thin 
C'^'^S emission. (Other methods of estimating the mass of 
these objects that do not assume virial balan ce, and make 
other assumptions about the emission, e.g. Plume et al.l 



IMueller et a] ■I2nna IRhirlev et a] .bnM^ve masses 
that are comparable to or better than the uncertainty in 
the mass estimate we adopt below.) They find a me- 
dian ratio Mh cn/-^0 = 6£hcn/ (K km s~^ pc^), and 
a mean of 11. iGao fc SolomonI (|2004aD perform radia- 
tive transfer calculations under a variety of assumptions 
about the temperature and distribution of the HCN- 
emitting gas, and arrive at conversion factors of 10 — 25. 
For normal spiral galaxies, we adopt a conversion of 
-^'^hcnZ-^O = 11 Lucn / (K k m s~"^ pc^), the mean value 
measured bv lWu et al.l l|2005j) . To be conservative we es- 
timate that the conversion is uncertain by a factor of 3, 
corresponding to a range that is larger than the range of 
6 — 25 spanned by the various estimates. 

Recent HCO+ observations suggest that chemical 
changes triggered by X-rays may en hance the HON abun- 
dance in galaxies with an AGN f Gracia-Carpio et al.l 
12006ft . This reduces the mass estimate for a given HON 
luminosity by a factor of ~ 2 in starbursts, most of which 
have at least a small AGN. For starbursts we therefore 
modify our light to mass conversion to Mhcn/-^© = 
5.5iHCN/ (K km s~ pc^), again with roughly a factor 
of 3 uncertainty. 

Combining this with the conversion from IR luminos- 
ity to star formation rate discussed in § 12.1.21 for both 
starburst galaxies and normal spirals gives M^/Mhcn ~ 
(30 Myr)~^. Adding the uncertainties of a factor of 2.5 
on the IR-HCN correlation, 3 on the HCN to mass con- 
version, and 1.3 on the IR to star formation rate conver- 
sion gives a combined uncertainty of a factor of 4.1. 

The free-fall time depends on the mean density, for 
which a rough guess based on the optical thickness of 
the HCN emission is uh ~ 6 x 10^ cm~^. The true den- 
sity of the HCN-emitting gas could be lower or higher 
depending on optical depth and beam-filling effects, as 
appears to be the case with CS(5-4) emission (see the 
discussion below in I2.4|l . However, this is unlikely be- 
cause the mass to luminosity ratio for HCN(l-O) inferred 
from radiative transfer calculations is very close to that 
inferred from correlation of HCN luminosity with virial 
masses for objects seen in the Galaxy. This suggests that 
the HCN emission for Galactic HCN clumps is close to 
beam-filling and is not extremely self-absorbed, so the 
density of the objects being observed is unlikely to be 
very different than uh ^ Q x 10^ cm~'^. If anything, the 
observations of W u et af] (|2005) suggest that the mean 
density may be slightly larger than uh ~ 6 x 10'* cm"'^ 
(also J. Wu, 2006, private communication). We adopt a 
mean density of nn ~ 6 x lO"* cm^"^ as a fiducial value, 
giving iff-HCN ~ 0.18 Myr, and estimate that the uncer- 
tainty is no more than a factor of 2, corresponding to a 
factor of 4 change in the density. 



Finally, since the extra-Galactic observations include 
all HCN(l-O) emission from the target galaxy, and all 
star-formation occurs at densities high enough to be 
traced by HCN emission, we set /hcn — 1- 

Using our fiducial estimates for M*/Mhcn, /hcnj and 
iff-HCN in CJ gives SFRff_HCN = 0.0058. The factor of 
4.1 uncertainty in our value of M*/Mhcn, combined with 
the factor of 2 uncertainty in our free-fall time, gives a 
complete uncertainty of a factor of 4.6. 

2.4. Star Formation in Gas Traced by CS 

The CS(5-4) line has a critical density of uh ~ 1.5 x 
10^ /r cm~^, and thus traces ga s at even hi gher densities 
than HCN(l-O). IPlume et all and IShirlev eTall 

l|2b03.) survey in CS emission lines a sample of regions 
selected from sources reported i n the Arcetri H2O maser 
catalog (jValdettaro et alJl200H) that are thought based 
on IRAS colors to be associated with star formation. 
From the mean CS(5-4) luminosity of their targets and 
the sky coverage fraction of the maser catalog, they es- 
timate that the total CS(5-4) luminosity of the Galaxy 
is at least Lcs ~ 20 L©. This is a lower limit only, be- 
cause the maser catalog may not be complete over the 
region of sky it covers, there is probably at least some 
Galactic CS(5-4 ) emission tha t is not correlated with wa- 
ter masers, and IShirlev et al.l l|2003r) only detect 75% of 
their targets in CS(5-4), while they detect 90% in all the 
transitions for which they look. Nonetheless, this lower 
limit on luminosity gives a lower limit on the amount of 
gas in the Galaxy that is sufficiently dense to produce 
CS(5-4) emission, which we can in turn use in equation 
iQJ to obtain an upper limit on SFRg at CS densities. 

First we must estimate the CS(5-4) "X" factor to con- 
vert the luminosity to a mass. Like HCN (1-0), star- 
forming clumps are usually optically thick in CS(5-4), so 
it is reasonable to expect that such an X factor might ex- 
ist. Figure|31shows the virial mass Mvir (determined from 
optically-t hin C^'^S emission) v ersus Lcs for the 57 ob- 
jects in the IShirlev et alJ l)2003f) sample. The correlation 
is weh fit by the line Mcs/Mq = 4.5 x W^Lcs/Lq over 
more than two orders of magnitude in mass. This implies 
a total Galactic mass in CS clumps of Mqs > 9 x 10^ 
Mq. The uncertainty of the fit, done in the logarithm of 
the data, is a factor of 3.3, and we adopt this as the un- 
certainty in our luminosity to mass correlation. However, 
as we discuss more below, the virial mass g ives the low- 
est mass of the three estimators tested bv iPhime et alJ 
(^97), by roughly a factor of 3, so the true mass is 
much more likely to be larger than our value rather than 
smaller. 

Note that this mass estimate strongly suggests that the 
true CS mass in the Milky Way is substantially larger. 
CS and HCN emission of ten come from o verlapping re- 
gions m the Milky Way l|Wu et all [20051) . and in such 
overlapping regions the area emitting CS is not vastly 
smaller than the area emitting HCN. However, the to- 
tal mass of CS gas that our estimate su ggests is smaller 
than the mass of HCN-emitting gas that IGao fc SolonionI 
Il2004albr see in galaxies like the Milky Way by 1 — 2 
orders of magnitude. The most likely explanation is in- 
completeness of the CS survey. 

To est imate the free-fall t ime, we must know the gas 
density. IPlume et al.l l)1997|) . based on radiative trans- 
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Fig. 3. — Virial mass Mvir versus luminosity in the CS(5-4) 
line L (-;r {asterisks with error bars) for objects in the Shirley ct all 
t2003l) sample. We also show the best-fit linear correlation [line). 
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Fig. 4. — Free-fall time an d density versus ma ss for the 57 
CS(5-4) clumps in the sample of IShirley et^ l2003l . 



fer modeling, find a mean density of = 1.6 x 10^ 
cm^"^, roughly the critical density. However, this aver- 
age is i ntensity- weighted, s o it r tiay overestini ate the true 
mean. iPlume et alT 1)19971) and 'Shirl ey et alJ 1)20031) find 
that the virial mass for objects in their survey is system- 
atically smaller than the mass estimated by assuming 
that all the gas is at the density inferred from the ra- 
diative transfer calculations. Based on the difference in 
mass estimates, they conclude that the filling factor of 
gas at densities of uh ~ 2 x 10^ cm~^ or higher is typ- 
ically ~ 0.3 - 0.5. Thus, the density of Uh = 1.6 x 10^ 
cm~^ probably represents an upper limit on the true den- 
sity. We therefore compute the mean density for the 57 
clumps in the Shirley et al. catalog from the virial mass 
and deconvolved half-peak radius (Shirley et al.'s Res) 
to compare to the result of the radiative transfer calcula- 
tions. We show our derived densities and the correspond- 
ing free-fall times in Figure 0] The mass-weighted har- 
monic mean free-fall time for the CS clumps (computed 
from equation^ is tff-cs = 0.10 Myr, and we adopt this 
as our fiducial value. The range 0.05 — 0.2 Myr covers the 
central 67% of the free-fall times, so we take the uncer- 
tainty to be a factor of 2. The density corresponding to 
our adopted mean free-fall time is nn = 1.8 x 10^ cm~^, 
and the range in free-fall times corresponds to a density 
range of a factor of 1.4. 

Finally, the estimated CS luminosity comes from 
sources associated with water masers, which arise only 



in regions of massive, clustered star formation. This im- 
plies a maximum of fcs — 0.8. Alternately, we would 
obtain the same result by assuming that all star-forming 
regions produce CS(5-4) emission and raising our esti- 
mate of Mqs by a factor of 1/0.8 to account for the frac- 
tion of CS(5-4) emission associated with water masers. 
In either case, our value for either fcs or Mqs is an upper 
limit because we do not know what fraction of clustered 
star formation is associated with water maser emission. 
Combining fcs, Afcs, and tg-cs in equation gives 
SFRff_cg < 0.27, with an uncertainty of a factor of 3.6. 
Again, it seems likely that the true value is significantly 
lower than this upper limit. 

2.5. Star Formation in the Orion Nebula Cluster 

We can add one more point for a specific object. The 
Orion Nebula Cluster (ONC) is the only star cluster that 
has been studied well enough so that it is possible to 
make reasonable estimates of the stellar mass, the for- 
mation time, and the properties of the progenitor gas 
system, allowing direct determination of SFRff. This 
is quite useful because it provides an estimate of SFRg 
that does not depend on conversions from luminosities 
to masses or star formation rates, and thus is subject to 
completely different systematic errors than the methods 
we have used thus far. 

The total stellar mass of t h e ONC is 4600 Mq 
l)Hillenbrand & Hartma'nnlll998D . iTan et all 1)20061) ana- 
lyze several lines of evidence to conclude that formation 
of the ONC took place over at least 3 Myr, and possi- 
bly longer - see Huff & Stabler (2006). The age estimate 
comes from a combination of pre-main sequence fitting, 
which is fairly reliable for stars ^ 2 Myr in age, and from 
the a ge of a dynamica l ejection event. We refer read- 
ers to iTan et al.l l)2006f) for the detailed evidence. The 
free-fall time in the current stellar system is tg = 0.5 
Myr, but it must have been smaller in the progenitor 
gas system, since some of the mass has been expelled 
by the Trapezium stars, and the cluster may be expand- 
ing. Uncertainties in these processes lead to a range of 
mass estimates for the progenitor, ranging from ~ 15000 
Mq if th e ONC today is ma rginally unbound and ex- 
panding l)Kroupa et alJl200l to 6700 Mg, if the ONC 
is cu rrently bound and non-expanding l)Huff fc StahleJ 
'2006"), although the latter would imply an extraordinar- 
ily high star formation efficiency. If we neglect the pos- 
sibility that the cluster has undergone significant expan- 
sion since expelling its gas, this implies that the free-fall 
time in the progenitor gas system was tg — 0.3 — 0.4 
Myr, corresponding to a density uh = 1 — 2 x 10** cm~^. 
It seems unlikely that star formation began when most 
of the gas in the ONC was spread over a much larger 
region at lower density, both because there is no correla- 
tion betwe en the age and spatia l distribution for stars in 
the ONC l)Huff fc Stahleill2006D and because there must 
have been a system sufficiently dense to lead to dynami- 
cal interactions and ejection of massive st ars 2.5 Myr ago 
l)Hoogerwerf et al.ll2f)0lt iTan et alJl20f)6l) . 

Combining the initial gas mass, final stellar mass, 
and initial free-fall time implies that in Orion SFRg = 
0.03 — 0.09. The true value is likely to be towards the 
low end of this range, since even if the ONC is bound it 
has almost certainly undergone some expansion from its 
original gaseous state, meaning that the free-fall time we 
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have adopted is probably too large. 

2.6. Summary of Observations 

We sunimarize our results by plotting the dimension- 
less star formation rate SFRg versus characteristic den- 
sity in Figure [S] In addition to our points for CS clumps, 
HCN clumps, IRDCs, and the ONC, we can add a 
point for GMCs as a whole, in which the typical den- 
sity a nd free- fall time are n — 100 cm~'^ and tff = 4.4 
Myr (|McKeei ri999). The total mass of GMCs in the 
Galaxy is roughly 10^ Mq ijBronfman et al.ll2000|) . so 
using the same argument as in S 12.21 gives SFRff w 0.013. 
IKrumholz fc McKeel ll2f)nl give a much more detailed 
calculation of this value, but for simplicity we adopt 
rough numbers here. Since the GMC population in the 
Milky Way is reasonably well characterized, the uncer- 
tainty in SFRff is dominated by uncertainty in the star 
formation rate, so it is about a factor of 2. The range 
of densities covered by star-forming GMCs is roughly a 
factor of 3, so we take this range and the range in SFRff 
to roughly define our uncertainty for GMCs. 

3. IMPLICATIONS AND COMPARISON TO THEORY 

The star formation rate per free-fall time is a key pre- 
diction of any model for a physical mechanism that regu- 
lates the star formation rate. The most important result 
summarized by Figure El which any successful theory 
must be able to explain, is that there is no evidence for 
a transition to rapid star formation out to densities of 
n ^ 10^ cm"-^. Even given the large uncertainties in the 
data, this conclusion is quite firm. If any of the objects 
we have analyzed were bound, detatched from the large- 
scale turbulent flow, and collapsing without significant 
impedance from feedback, we would expect SFRff ~ 1, 
while even the top of the uncertainty range for the data 
point with the highest value of SFRff is below 0.1. Only 
the CS data is consistent with a transition for rapid star 
formation, and this is an upper limit for which there is 
strong evidence to suggest that the true value is signifi- 
cantly lower. 

A second conclusion, which is much more tentative, is 
that we do not see any significant change in SFRff with 
density at all. Given the uncertainties, obviously we are 
insensitive to changes in SFRff smaller than roughly an 
order of magnitude, and the somewhat higher point for 
the ONC suggests that it is possible there is some modest 
rise in SFRff in star clusters. Any such rise with density, 
though, must be fairly small over the range of density we 
have explored. In summary, the data are consistent with 
SFRff of a few percent regardless of density, and require 
SFRff < 0.1 to densities of at least 6 x lO'' cm-^. The 
characteristic scale at which one transitions from slow to 
fast star formation must lie at even higher densities. Here 
we investigate how well various models of star formation 
explain this result, and discuss a few of its implications. 

3.1. Unbound GMC / Collapsing Cluster Models 

One proposed explanation for the low star formation 
rate in mol ecular c louds is that GMCs arc not bound 
by gravitv llMac Lo w & Klcsscn 2004; Clark & BonncU 
l2f)0l IClark et a,lJ 1200,^ IDobbs et a,].l l200(ii) . In this 
model, GMCs are transient over-densities created by tur- 
bulence or spiral shocks in the atomic ISM. Since they are 



confined by ram pressure, they re-expand and disperse 
in roughly one dynamical time. The star formation rate 
appears low because only a small fraction of the mass 
is gravitationally bound and can collapse, but that mass 
which can collapse will form stars rapidly. Feedback is 
not effective at slowing collapse or establishing an equi- 
librium, and only stops the star formati on process at its 
end by rapidly unbinding the cluster (jBonnell fc Batd 
120061) . Star formation in protoclusters is rapid but is 
halted after a short period by feedback. Thus, the star 
formation rate is high, even if the duration is short, and 
this should be detectable as a rise in SFRff at the typical 
density of protoclusters. 

Simulations based on these premises allow us to com- 
pare these models to the observational data. Before dis- 
cussing individual simulations, we note that we cannot 
analyze the simulation data as thoroughly as we have an- 
alyzed the observational data. Ideally, one would take a 
snapshot of a simulation at various points in time, com- 
pute the mass M(> p) that is denser than p, make a plot 
of A/(> p)/tf[{p), and compare that to the instantaneous 
star formation rate in the simulation. This is effectively 
what we have done with the observational data. How- 
ever, simulation papers generally do not report M{> p). 
They only give a single density, usually the initial (uni- 
form) density, and report a total mass of stars formed 
over the course of the entire simulation. Thus we per- 
form our calculations using M(> p) for the initial density 
at the initial time, and compare the mean star formation 
rate over the entire simulation. This should give a rough 
average value that can stand in for a more detailed cal- 
culation of M{> p)/ts{p) as a function of density and 
time. However, we encourage authors of future simu- 
lations of GMC evolution or star cluster formation to 
compute M{> p)/tff{p) for at least a few time slices to 
f acilitate co mparison to observations. 

IClark et' al. ( 2005) estimate from their simulations that 
unbound GMCs would convert 5 — 10% of their mass into 
stars in 2 — 3 free-fal l times , giving SFRff ~ 0.03. In con- 
trast, ^^nnsJOOD 1)20031) simulate a marginally-bound 
1000 Mq clump with an initial radius of 0.5 pc, giving 
uh — 5.5 X 10^ cm~'^, roughly the properties of an HCN 
clump. There is no feedback, and the turbulence decays 
freely. They find that, after 2.6 initial free-fall times, 58% 
of the mass has gone into stars, giving SFRff = 0.22. This 
is probably an underestimate of the true value of SFRff in 
this scenario, since the simulation starts with a uniform 
density and thus star formation does not really begin for 
roughly 1 free-fall time. If we instead measured from the 
time when the first star formed to the time when the last 
one formed, we would find an even larger value of SFRff. 

This is consistent with other simulatio ns in which the 
turbulence is weak. For example, .Klessen fc BurkerO 
(2000) model formation of a star cluster by simulating a 
periodic box in which the gravitational potential energy 
greatly exceeds the kinetic energy, and find that their 
simulation converts 60% of the initial mass into stars in 
1.8 free-fall times, giving SFRff = 0.33. The simulation is 
scale-free, but Klessen & Burkert suggest that the model 
should be reasonab le for a re gion wi th a typ ical density 
of uh ^ 10^ cm~'^. iVaza^iez' Scmadc ni et al.l (^003) per- 
form a parameter study of regions with turbulent driv- 
ing of various strengths and find that when the turbulent 
driving is weak so that the system is unstable to global 
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Fig. 5. — Estimates of the dimensionless star formation rate SFRff in objects of characteristic density nu- We show observational 
estimates {filled squares with error bars), an observational upper limit [filled square with arrow), simulations with weak or decaying 
turbulence (open diamonds), simulations with strong, driven turbulence {filled diamonds), and the prediction of an analytic model {hatched 
region). For the observational points, we label each point by the tracer used to estimate that point. We label t he simulation points and 
the analytic model region by the reference for the simulation or model. The simulations of'Klessen et al." f 2000*1. 'Vazguez-Semade ni et all 
©03), and Li et al. (2004)' are dimensionless, so the values of njj for them are arbitrary. We have adopted the authors' suggested scale 
value of uh where one is given, or an arbitrary value otherwise. For the simulations of Vazgue z-Semadeni et al .. we show a weakly-driven 
model (M2K2, meaning Mach number 2, driving wavelength of 1/2 the box size), and a strongly-driven one (M10K8), and wc derive values 
of SFRfj for these models via the method described in .Krumholz &: McKed I2005i) . The region indicated for the tKrumholz fc McKea model 
corresponds to the parameter range Ovir = 1 — 2, Ai = 20 — 40. 



collapse, SFRff ~ 0.3 is a typical value. (As we discuss 
below in § 13.41 iVazauez-Semadeni et al.ll2003l find that if 
the turbulent driving is strong, the simulations produce 
lower values of SFRff that agree much better with the 
observations.) 

We place po ii its fr om the simulation s of 
Klessen fc BurkertI l|2000|) . IBonnell et all pOO^ . and 
c!ar^^t'ldTH2005[l . and from one of the weakly-driven 
turbulence simulations of ' Vazauez-Semadeni et aP 
lf2003) . on Figure El This indicates a clear observational 
problem for these models. They produce roughly the 
correct star formation rate at the GMC scale, but are 
inconsistent with the data on the IRDC and HCN scale 
at the order of magnitude level. They are consistent 
with the CS data only if the value of SFRff in CS clumps 
is near its upper limit. We conclude that models where 
star formation occurs in unbound GMCs but freely 
collapsing dense clumps, or where the turbulence in 
those clumps is driven but at very sub-virial levels, are 
ruled out by the data. Star formation in dense clumps 
cannot occur through global collapse. 

We emphasize again that one cannot avoid this prob- 
lem by hypothesizing a feedback process that rapidly de- 
stroys clumps once they have turned a relatively small 



fraction of their mass into stars, but retaining the pic- 
ture of star formation occuring in clumps that are un- 
dergoing global collapse. Even if feedback destroyed col- 
lapsing clumps, while they existed they would still be 
forming stars at a rate much larger than SFRff of a 
few percent, which is inconsistent with the data. Only 
if destruction by feedback prevents global collapse from 
starting can it produce the observed value of SFRff. In 
this case, though, estimates that cluster-forming clumps 
turn ~ 30% o f their gas into stars (Krouoa et al. 200^ 
ILada & Lad^ l2003^ require that global collapse be held 
off for much longer than a free-fall time. We discuss the 
implications of this in more detail in § 13.51 

One final caveat worth mentioning in our interpreta- 
tion of these simulations is that once a significant fraction 
of the particles in an SPH simulation have been accreted 
by sink particles, the effective resolution of that code, 
and its ability to follow the hydrodynamics, may be sig- 
nificantly degraded. It is unclear how large this effect is, 
and whether it would tend to increase or decrease SFRff. 
However, since the effect is only a problem late in the 
simulations, it seems unlikely that it could shift SFRff 
by more than a few percent. 
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3.2. Galactic- Scale Gravitational Instability Models 

The criticism of unbound GMC models applies in a 
much weaker form to models in which the star formation 
rate is determined by large-scale global gravitational in- 
stability in a galactic disk (e.g. iLi et alJl20'05albl . 120061: 
iTasker fc Brvanll2006j) . These models postulate that the 
star formation rate in galaxies is determined by the rate 
at which gravitational instability in the galactic disk 
drives gas to high densities, and that the strength of 
this gravitational instability determines the star forma- 
tion rate. This is probably at least partly true, in that the 
star formation rate is set partly by the supply of molec- 
ular gas, which is in turn controlled by gravitational in- 
stability. However, roughly one third of the mass of the 
Milky Way's ISM is in molecular clouds. Clearly the 
rate-limiting step in star formation is not at the stage of 
converting atomic ISM to molecular gas, but somewhere 
at higher density. For this reason, the gravitational in- 
stability models provide only a partial description of the 
star formation process. 

In practice, this limitation is reflected in the star for- 
mation prescriptions that simulations of galactic-scale 
gravitational instability adopt. These simulations can- 
not resolve structures of molecular cloud size or smaller. 
Instead, when such structures form the simulations use 
a recipe to determine the amount of gas that turns into 
stars. These prescriptions only reproduce the correct star 
formation rate if they correspond to rates SFRg 1 in 
the dense gas. 

For example, iTasker fc BrvanI l)2006j) assume for their 
highest resolution simulations that in cells with densi- 
ties above 10^ cm^'^, 50% of the gas is converted into 
stars per free-fall time of that gas. This corresponds to 
SFRff = 0.5. However, examining their Figure 8 shows 
that this leads to a star formation ra te that is an order 
of magnitude larger than the observed lKennicutil l)1998bf) 
law without feedback, and a half an order of magnitude 
too large if they include feedback. They also find that 
even reducing their value of SFRg by a factor of 10 is not 
by itself sufhcient to reproduce the Kennicutt law. This 
is not surprising based on our results that SFRg is at 
most a few percent in gas at densities of 10"^ cni^^. Re- 
producing the observed star formation law would require 
a value for SFRff of th is order. 

In the simulations of iLi et all ('2005a'b', 2006), the au- 
thors are able to reproduce the iKennicutt L1.9.98b,) law 
by adopting a prescription for star formation whereby 
collapsed mass is assumed to form stars with a 30% "ef- 
ficiency", in the sense that they count 30% of the gas 
that is accreted onto sink particles in the simulation as 
stars. The remaining 70% of the mass stays in the sink 
particles, but is assumed not to form stars. There is no 
simple way to convert this recipe into a value of SFRff, 
since it specifies that the star formation rate in any par- 
ticle is 30% of the particle's current accretion rate, and 
we do not have available the full time history of accre- 
tion for each particle. Even if we did have this infor- 
mation, interpretation would be diffuclt because the gas 
assumed not to form stars remains locked up in sink par- 
ticles indefinitely, rather than being ejected back into 
the ISM by feedback as happens in real star clusters. 
This greatly reduces the supply of gas available to form 
stars. Nonetheless, the fact that these simulations need 



to adopt a subgrid model that significantly reduces the 
collapsed m ass counted as stars in order to reproduce the 
IKennicutt law is consistent with our point that one can 
only explain the observed rates of star formation if there 
is some mechanism that inhibits star formation even in 
the densest gas that the simulations can resolve. 

The reason that these simulations reproduce the ob- 
served star formation rate only if they have SFRff ^ 1 
is clear in light of Figure [S] if the gas in IRDCs or HCN 
gas were collapsing on its free-fall time, the star forma- 
tion rate in the Milky Way would exceed the observed 
value by two orders of magnitude. To obtain the cor- 
rect overall star formation rate, one must adopt a star 
formation rate SFRff <C 1 even in the densest gas that 
the simulations are currently capable of resolving. While 
this is an entirely reasonable approach in simulations, a 
theoretical understanding of the star formation rate re- 
quires an explanation why SFRff <C 1 even in gas that is 
already collapsed to densities above the critical density 
of HCN (1-0), five orders of magnitude higher than the 
mean in the Galactic ISM. This point does not mean 
that global gravitational instability is unimportant in 
regulating star formation, simply that it cannot be the 
sole agent. Thus, models such as these are complemen- 
tary to models of SFRff in dense gas such as those of 
iKrumholz fc McKee (200l). 

3.3. Magnetic Regulation Models 

Another possible explanation for low st ar forma- 
tion rates is st rong m agnetic fields (e.g. iShu et alJ 
Tfl^ McKeA mm lTa,ssis Mouschovia,.cJ 12001 
Nakamura fc: Lil |2005|) . If star- forming clouds are 



magnetically subcritical, which is controversial o n both 
observational and theoreti cal grounds ( McKee et al 
T993[ [Crutcheri 1999: Bour keet al.ll200ll iPadoan et a" 
2003 



Heiles fc Crutcher. .2005j) . then they cannot col 



lapse before ambipolar diffusion allows the magnetic 
field to slip out of the gas, and star formation proceeds 
on the ambipolar diffusion time scale instead of the 
free-fall time scale. For a uniform gas, the time scale 
required for ambipola r diffusion to decouple the gas and 
the field is l|Shulll992t) 



^AD 



—iCp 



1/2 



V 



(6) 



A 



where L and p are the characteristic size scale and den- 
sity of the cloud, va is the Alfven speed, 7 ss 3.5 x 10^'^ 
cm-^ g-i g-i is the ion-neutral drag coefficient, and C — 
3 X 10""'^^ cm~^/^ g^/^ is the cosmic ray ionization con- 
stant. (In high column density environments like IRDCs, 
interstellar UV photons cannot penetrate and ionization 
is dominated by cosmic rays.) If we assume that this 
describes the star formation time scala^ then the star for 



mation rate should obey roughly ijCiolek fc MouschoviasI 
1I993,) 

SFRff ^ ^H^, (7) 

where ec ~ 0.5 is the fraction of the mass in a prestellar 
core that reaches the star rathe r than being ejected by 
outflows (|Matzner fc Mc Kcc 2000). 

To determine what this implies for SFRff, we need 
to know how the magnetic field behaves in objects of 
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varying densities. While there are very few direct ob- 
servations of magnetic fields for extremely dense ob- 
jects of the sort we are considering, observations do 
show that magnetic fields obey the correlation B = 
{87r)-'^/^crp^/^/i~^c}^^, where a is the velocity dispersion 
of the region, /z is the ratio of the object's mass-to-flux 
ratio to the critical value {2itG^^'^)~^ , c\ is a constant 
of order unity the depends on the cloud's internal den- 
sity distribution, and [I'^cl"^ « 0.8 ijBasul 1200(1) . This 
correlation is what one would expect if, in one direction, 
cloud self-gravity were balanced by magnetic plus tur- 
bulent pressure. If we use this correlation in Q, and 
re-write the relation in terms of the virial parameter 
ttvir = 5(T^L/(G'M), we find that all the dependence on 
dimensional quantities drops out and we are left with 



SFRff - O.Olttv 



(8) 



The uncertainty of this calculation is probably more 
than an order of magnitude, so we should not pay par- 
ticular attention to the coefhcient, and we will not at- 
tempt to place points for magnetic regulation models on 
Figure Major contributors to the uncertainty are am- 
biguities in the definition of the length and mass scales L 
and Af , the lack of three-dimensional numerical simula- 
tions to determine how well equation (Tjl holds on scales 
larger than a core collapsing to form a single star system, 
and the possibility that the cosmic ray ionization rate 
may vary substantially between galaxies. Perhaps most 
significantly, equation lO assumes that ambipolar dif- 
fusion in a turbulent medium is not substantially faster 
than in a quiescent one, an assumption that may well 
fail fe.g. iHeitsch et alJl20?)l INakamura & Lill2005|) . If 
turbulence shortens the ambipolar diffusion time scale 
so that ^AD iff, as some simulations suggest it might, 
then magnetic regulation models would produce a value 
of SFRff that is too high. Whether this occurs or not is 
a question that will have to be addressed through three- 
dimensional simulations with self-gravity and non-ideal 
MHD. 

Given these uncertainties, we cannot really say 
whether the rate of star formation in magnetically sub- 
critical clouds subject to ambipolar diffusion is quanti- 
tatively consistent with the observations. However, the 
lack of dependence of SFRff on any properties but the 
virial parameter implies that we expect the magnetic 
SFRff to be roughly the same in all virialized objects, a 
prediction we can compare to observations. For our ob- 
served objects, we know that GMCs with masses ^ 10^ 
MfT), HCN clumps, and CS clumps arc roughly virialized 
l|Phime et al. 1997; Hcvcr ct al. 2001; Wu ct al. 2005), so 
avir ~ 1 for them. For IRDCs we lack kinematic informa- 
tion from optically thin molecular emission, and therefore 
we cannot directly determine the velocity dispersion and 
the virial parameter. If we assume that these objects are 
virialized, then the magnetic regulation model is broadly 
consistent with observations that SFRff is roughly con- 
stant. We cannot make a stronger statement than this 
because the magnetic regulation model cannot currently 
make more specific predictions. Note that it is not criti- 
cal for this purpose that the objects in question truly be 
in exact virial equilibrium, since they almost certainly 
are not. We are using the virial parameter only as a way 
of parameterizing the strength of the turbulence relative 



to the strength of gravity, so all that we require is that 
the objects have virial parameters avir ~ 1, i.e. that 
they not be either completely collapsing or completely 
non-self-gravitating. 

3.4. Virialized Turbulence Models 

A fourth idea to explain the low star forma- 
tion rate is that turbulence inhibits collapse. 
This id ea has a long hist o ry, an d we refer read- 
ers to IMac Low fc KlessenI (2(M)_ for a thorough 
review. Recently. iKrumholz fc McKed l)2005j) syn- 
thesized simulations and analysis of turbulence- 
inhibited star formation by 'Kl essen et alJ l)2000t) and 
fVazqucz-Scmadeni et al. (2003) to derive an esti- 
mate that turbulence produces a star formation rate 
SFRff « 0.014(avir/l-3)-° '^*(7U/100)'°-32, where Ovir 
and A4 are the virial parameter and Mach number of 
a star-forming gas cloud. The estimate is based on 
a derivation of the fraction of mass that is unstable 
to collapse in a medium that has the density and 
velocity structure common to supersonic isothermal 
turbulence, and is calibrated against simulations. Its 
uncertainty is probably a factor of a few, stemming from 
uncertainty in the effects of magnetic fields and from 
the uncertain approximation ec = 0.5. For GMCs, HCN 



clumps, and CS clu mps, typical values ar e a 



1 



and - 20 - 40 ISolomon et abllTflSTt iMi^lTgga 
IPlume et aljlT997t IShirlev et al.ll2003t IWu et all 120051) . 

As with the magentic models, exact virial equilibrium 
is not required, simply that the kinetic and potential 
energies be comparable. We plot the predicted value of 
SFRff for this range of parameters in Figure \E\ Given 
the uncertainties in both the observations and the 
theoretical calculation, there is reasonable agreement. 

Simulations of star formation in turbulent gas for 
which the turbulent energy is roughly equal to the gravi- 
tational potential energy, either because it is driven arti- 
ficially or because the simulations include feedback from 
protostellar outflows, also give values of SFRff that are 
broadly consistent with the data. We show three ex- 
amples in Figure Li ct al. (2004) simulate a periodic 
box in which the kinetic and turbulent energies are ap- 
proximately equal and the turbulence is driven to keep 
the level of turbulence roughly constant. They find a 
star formation rate of SFRff = 0.079 at a resolution of 
512"^ cells. (Since the simulation uses periodic bound- 
ary conditions, the true gravitational potential energy 
is not well-defined, so we cannot state that the virial- 
ization is more than approximate.) The simulations are 
dimensionless, but they suggest scaling to a typical den- 
sity of 2 X 10^ cm""^. Similarly^n sc ale-free sim ulations 
with stron g turbulent driving. IVa zauez- Semadeni et alJ 
1)2003(1 andlVazauez-Semadeni etaL '( 2005*) find SFRff - 
0.01. ILi fc Nakamiiral mO(t\ fi~SFR,ff = 0.067 in a 
simulation with a central density of 5 x 10* cm"'^ that 
includes protostellar outflows from forming stars, which 
keep the kinetic and gravitational potential energies ap- 
proximately equal. A number of other simulations in 
which the turbulence is roughly in virial balance with 
the gravitational energy produce star formation rates in 
reasonable agreement with observations over the period 
in the simulations for which that balance is maintained. 
If the turbulence is sufficientl y strong, simulations can 
achieve SFRff 0.01 (e.g. IVazauez-Semadeni et alJ 



12 



200: 
20Q5D. 



iTillev fc Pudritd 120041: iVazauez-Semadeiii et all 



We note that SFRff by itself does not appear to pro- 
vide a way of distinguishing the length scale at which 
turbulence is driven, since simulations in which tur- 
bulence is driven primarily at large scales and ones 
in which it is driven primarily at small scales both 
seem capable of producing low values of SF Rff. Some 
authors (e.g. Vazaucz-Scmadcn i et al.ll2003|) emphasize 
the length scale of the driving as critical to determin- 
ing whether the turbulenc e can keep SFRg low, and 
IKrumholz fc McKeel l|2flfl5f> show that the characteris- 
tic distinguishing driving that gives SFRg ~ 1 from 
driving that produces SFRff ^ 1 is how the Jeans 
length in the gas compares to the sonic length of the 
turbulence, defined as the length scale for which the 
turbulent velocity dispersion equals the sound speed. 
Other diagnost ics, such as the IMF and clustering statis- 
tics (e.g. Klcs sen fc Burker3 120011 'Schmcia & Klcsse^ 
2004"), the linewidth-size relation of the gas (e.g. 
Ossenl«3pLfc_Jyiac .Low 2002,), and the line luminosity of 



molecular clouds ljHever fc Brunti2004|) suggest that tur- 
bulence in molecula r clouds is primarily d riven on large 
scales. In this case IKrumholz fc McKeel show that the 
condition distinguishing cases with SFRff ^ 1 from those 
with SFRff <C 1 reduces to the statement that clouds 
with ttvir ^ 1 have low values of SFRff, while those with 
Qfvir ^ 1 have SFRff ~ 1. Thus, the combination of 
observations of SFRff and those showing that most tur- 
bulent energy is on large scales imply that turbulence 
can only reproduce the observed values of SFRff if it is 
driven strongly enough to produce a^h- ^ 1. 

3.5. The Formation Timescale of Clusters 

Another important implication of this work is for the 
formation timescale of star clusters. The star formation 
efficiency of clusters, defined as the fraction of the initial 
gas mass that forms stars over the lifeti me of the gas 
clump, is thought to be ~ 20% — 50% l)KrouDa et al.l 
|2001; Lada fc Lada1 l2003|) . Since the observational data 
show that SFRff is at most a few percent, this implies 
that the cluster formation process must take at least ~ 
lOtff (~ 5 crossing times). 

The recent observation by iWu et aP l)2005t) that in- 
dividual HCN clumps in the Galaxy with IR luminosi- 
ties Liji > 10^'^ Lq lie on the same Ljr — Xhcn corre- 
lation as entire galaxies also provides indirect evidence 
that the cluster formation time scale is comparatively 
long. Wu et al. suggest that luminous HCN clumps fall 
on the galactic IR-HCN correlation because they con- 
tain enough stars to sample the IMF fully. Our results 
showing a large scatter in the light output for small clus- 
ters support this conjecture. However, our results on the 
age-dependence of the luminosity suggest that a mass 
large enough to sample the IMF is not sufhcient by it- 
self to place a cluster on the extragalactic IR-HCN cor- 
relation. Since star formation in the galaxies surveyed 
by IGao fc SolomonI lj2004albfl has been ongoing for }^ 3 
Myr, their IR luminosities trace the star formation rate, 
and the Lm — Lhcn correlation is therefore a statement 
of the star formation rate per unit ma ss in HCN gas. I f 
the Galactic HCN clumps observed bv lWu et all 1)20051) 
typically survived for only ~ 1 free-fall time, ~ 0.2 Myr, 



then as Figure^shows their luminosity per unit star for- 
mation rate would be quite different than that of an older 
galactic population, and they would not lie on the galac- 
tic correlation. One would only expect Galactic HCN 
chimps to follow the extragalactic Lm — ^hcn correla- 
tion if they are massive enough to sample the IMF well 
and if they are ^ 1 Myr old, which means that they must 
have existed for at least ~ 5ts, and probably closer to 
lots. 

These two lines of evidence provide additional sup- 
port for the argument that star clus ters form in ncar- 
equihbrium gas clumps presented bv iTan e'Fa l. ( 2006). 
Tan et al. discuss the implications of this finding in 
more detail, but we note here two of the most signifi- 
cant. First, formation times of ~ lOiff or longer appear 
inconsistent with models in which turbulence in proto- 
cluster gas is freely decaying. Simulation s show that if 
turbulence is allowed to decay freely (e.g. 'BonncU et ajj 
or never contained energ y comparabl e to the grav- 
itational potential energy (e.g. IVazouez-Semadeni et all 
.2003,) . the star formation process generally ends in only 
about 2ts, with ^ 50% of the mass in stars. Only if 
the turbulence is continually driven can there be enough 
gas left up to lOtff for the star forma tion process to con- 
tinue (e.g.l Vazaucz-Sc madeni et aLll2003HLi et al.ll2004l: 
iLi fc Nakamura 2006) . The long cluster formation time 
scales suggested by the observations therefore imply that 
either the simulations of decaying turbulence are incor- 
rect or (more likely) the turbulence is con tinuallv driven. 
This is cons istent with on observations bv lWilliams et af] 
p003i) and lOuillen et al.l (12005,) , and with simulations 
of protostellar feedback bv lLi fc Nakamura ( 2006). We 
note again that the observations of SFRff alone put no 
constraints on how or at what scale the turbulence is 
driven. It can be from internal sources of feedback on 
small scales, or from a turbulence cascade in which most 
of the power is on the size scale of an entire GMC. 

A second implication is for the mechanism responsi- 
ble for the stellar initial mass function, ^rumholz et a3 
(12005*) have recently emphasized that the competitive ac- 
cretion mechanism (Bon nell et aL.2006n, and references 
therein) can only operate in the context of a strongly 
sub-virial gas clump that is undergoing global collapse 
and converting order unity of its gas m ass into stars in 
a free-fall time. iBonnell fc Batd ()2006 ^ concur that the 
question of whether competitive accretion occurs or not 
depends on the extent to which gas and stellar veloci- 
ties in a forming star cluster remain tightly coupled for 
the majority of the formation process, which in turn de- 
pends on whether feedback is effective at driving turbu- 
lence, which would change gas velocities but not stellar 
velocities. If turbulent driving keeps clusters virialized 
for many free-fall times and limits SFRff to ^0.1, then 
competitive accretion cannot occur within those clusters. 
The observations we discuss here provide evidence for low 
values of SFRff and long cluster formation times, and 
therefore provide strong evidence against the possiblity 
that competitive accretion determines the IMF. 

4. CONCLUSIONS 

We present observational evidence for two surprising 
conclusions, one shown quite strongly and the other more 
tentatively. First, and very clearly, star formation in 
dense gas is slow. The time required to convert all the 
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gas into stars, the depletion time, is longer than the 
free-fall time by at least and order of magnitude, and 
probably closer to two orders of magnitude. There is no 
evidence that gas at densities up to uh ~ 10^ cm~'^ re- 
sides primarily in objects that are undergoing collapse 
and rapid star formation. Second, and more tentatively, 
the ratio of free-fall time to the depletion time is inde- 
pendent of the characteristic density of the star-forming 
object in question. The data at present are uncertain, 
and are insensitive to changes of less than roughly an 
order of magnitude. Nonetheless, it is interesting that 
this ratio varies by no more than an order of magnitude 
from the typical density of GMCs to the typical density 
of HCN(l-0)-emitting gas, a range of nearly 3 decades in 
density. 

These observations are a strong constraint for theo- 
ries of star formation, and seem difficult to explain in 
the context of models in which there is a transition from 
slow, unbound star formation to rapid, bound star for- 
mation somewhere between the GMC scale and the pro- 
tocluster scale. Slow star formation in cluster-forming 
gas also implies that clusters requir e many free - fall ti mes 
to assemble, as recently argued by iTan et iahl 1)200(11) on 
other grounds. Models in which star formation takes 
place in virialized objects and is inhibited by strong mag- 
netic fields are qualitatively consistent with the data, and 
models in which star formation is inhibited by turbulence 
are both qualitatively and quantitatively consistent with 
observations. 

In the future it would be extremely useful to improve 
the data on which Figure O is based. One way to do this 
would be to perform detailed studies of other young clus- 
ters and obtain data comparable to that for the ONC. 
This would provide a method of estimating SFRff that 



is independent of luminosity conversions and does not 
suffer from concerns about the completeness of galactic 
surveys. Another improvement in the data could come 
from an unbiased survey of CS(5-4) emission in the Milky 
Way or in another galaxy, which would allow us to re- 
place the upper limit on SFRff we derive here with an 
actual estimate. This would be particularly valuable be- 
cause the upper limit on SFRff from CS(5-4) is well above 
the estimate from HCN(l-O), even though the densities 
are not very different. Unbiased CS observations could 
likely bring down this point. 

Another improvement would be extend the data to 
higher densities. To accomplish this will require obser- 
vations either of external galaxies or relatively complete 
surveys of the Milky Way in molecular transitions that 
trace densities ^ 10® cm~'^. Determining masses from 
the luminosities in these transitions will probably require 
high resolution follow-up observations of Galactic sources 
in optically thin isotopomcrs so that the luminosity may 
be correlated against a virial mass. While this is a signif- 
icant observational challenge, such surveys might make it 
possible to identify a scale at which star formation transi- 
tions from slow to fast, a crucial datum in understanding 
the star formation process. 
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